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The  energy  averaging  of  the  3o Itzmann- Vlasov  icine  + !o  ^qu.tion 
is  performed  rur  small  discharge  currents  in  a one-dimensional 
approximation  on  the  assumption  tha  \ the  velocity  distribution 
of  the  electrons  is  near  to  the  Maxwell  distribution. 


1 . INTRODUCTION 

One  of  the  starting  equations  of  the  ionizing  waves*  theory 
developed  at  the  Physics  Institute  CSAV  in  Prague  is  a linearized 
energy  transfer  equation  (see  /l / ) 


(1) 


I8  + = bt« , 

dz 


where  "z"  is  the  discharge  axis  (its  positive  direction  here  is 

from  the  cathode  to  the  anode  - for  this  reason  in  the  equation 

^ thermal 

( 1 ) the  signs  are  changed),^  and  "e"  are  electron 

deviations  expressed  in  energy  units  and  longitudinal  electrical 

field  intensities  ' deviations  from  a state  of  equilibrium. 

This  equation  was  obtained  by  linearizing  the  Granovs'cy  /2/ 

equation  for  the  energy  balance  of  electron  collisions. 

In  this  wor'c  the  energy  transfer  equation  is  exactly  deduced 

from  the  Boltzmann-Vlasov  tlnetic  equation  for  electrons 


(2) 


cL 

dt 


for  a case  of  quasl-Maxwellian  electron  velocities  distribution 
for  small  discharge  currents  in  a uni-d  imensional  approximation 
(not  applied  to  non-homogeneous  elements  in  the  radial  direction). 
In  the  equation  (2)  P is  a partition  function  of  the  electrons 
velocities,  "t"  is  time,  "c"  is  instantaneous  velocity,  Vr,  Vc 
are  operators  in  the  coordinates  and  velocities'  sub-space  of 
the  phase  space,  Fe  is  the  force  acting  upon  the  electrons,  mp 


is  the  mans  of  electron. 
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the  collision  member  which 
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expresses  liar  the  interaction  of  electron?!  among  themneLves  and 
with  all  the  other  particlea. 


2.  DEDUCTION  OF  THF  EQUATION  FOR  THE  TRANSFER 

ENERGY  OF  ELECTRONS. 

i n 

Under  laboratory  condltiona  low  preaaure  plasma  Let  ua 
aaaumr  that  the  partition  function  fp  appears  in  the  following 

manner  (aa  in  reference  /3/) 

(3)  /,  *=  K,  exp  {-y,»!,[(c,  - u.y  + c\  + cj]}  ; 


V _ I /j  P T~ 

e , ”'c"  ia  the  Boltzmann  conatant  and  Tp  ia  the  electron 

temperature,  cj^  ia  the  i-th  eLement  of  inatantaneoua  velocity  and 
u^  ia  the  "carrying  away"  velocity  of  the  electrona  (in  the 
direction  of  the  "z"  axis).  The  conatant  ia  introduced  as  a 


standard  for  the  concentration  N 


. f f.dC. 

J(c) 


Using  the  method  of  reference  /j/  it  ia  possible  to  demonstrate 


that 

(4) 


In  order  to  obtain  the  electrons  energy  transfer  equations,  we 
have  to  multipLy  equation  (2)  by  a quadruple  of  the  inatantaneoua 
veLocity  "c"  and  integrate  throughout  the  velocities'  sub-space  in 
the  phase  space  (see  also  /3/)» 

(5)  L'%‘c+L«'A/.*e+isfcty«- 

The  first  two  integrals  on  the  left  side  of  this  equation  may  be 
calculated  sLmply  by  expressing  the  instantaneous  velocities  in 
spheric  coordinated,  carrying  out  the  suitable  substitutions,  deve  . 
loping  in  series  and  by  Limiting  (with  a precision  of  3 significant 
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figures)  ourseLves  to  its  first  member  (see  addendum). 

ihe 

The  integration  method  chosen  explains  bain  neglect  of  the  kine- 
tic energy  of  the  "carrying  away"  veLocity  of  the  electrons  as 

e 

compared  to  the  thermal  energy.  If  we  denote 


= 0, 


we  shall  obtain 


• ' v.  n i kb  v w WU  k » i 

(7)  f ci%dC  = j 2~K0, 

' ' ' J(r)  Cl  dim, 

(R)  f c\cV,)f,dC~  yi  l-S,Ou,. 

J(r)  OZ  3 m. 

When  calculating  the  third  integral  we  have  to  maice  an 
assumption  about  the  acting  forcesi  in  the  case  of  the  usual 
conditions  in  the  positive  column  glowing  discharge  Fp  in  the 
T.orentz  force  given  by  the  equation 

(9  ) F,  = -q0E  - --  [c  x H]  , 

fo 

where  -q  is  the  eLectron  charge,  "F"  and  "H"  are  intensities  of 
the  electrical  and  magnetic  fields,  C0  is  the  velocity  of  light 
in  vacuum.  By  a procedure  similar  to  the  calculation  of  the  above 
mentioned  integrals  we  may  easily  obtain  the  following  (see  the 


addendum ) s 


f A - V«U  dC  = - — <JoN>«£ 


The  integral  on  the  right  side  expresses  the  cooLing  or  heating 
of  electrons  by  the  action  of  collisions.  T,ater  we  will  provide 
tnis  equation  with  an  approximation  expression  and  for  this  purpose 
we  will  establish  which  processes  we  shall  taice  into  account 
and  which  we  shall  neglect.  In  the  meantime  let  us  point  outi 


f - 

J„.  l6'J 


If  we  substitute  the  expressions  (7),  (8),  ( 1 0 ) art  (n)  in  (5) 
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and  solving  for  mp/2  we  obtain  the  following  equation 

( 12  } ~ (A',©)  + * - - {N.Ou.)  = q0.V.w,£  + -4  • 

Cl  3 07 

By  further  arrangements  and  with  the  aid  of  the  electrons?  conti- 


nuity equation 


dN,  . v 5 N, 
— * + div  (N,uc)  = 

6/ 


(where  ^ 'v<  expresses  the  creation  and  extintion  of  electrons 

St." 

by  ionization  and  recombination  respectively)  , we  obtain  the  final 
energy  transfer  equations 

, iv  dO  cO  1 dN  O 1 RV  1 

14)  — + m,  — - = - K — - io  — - JL  0 /l . 

fa  fa  N.  dz  dz  N,  bt  N, 

Its*  right  side  we  interpret  as  followsi  The  first  member  expresses 
heating  by  Joule  tyoe  heat,  the  second  expresses  cooling  (or  heating) 
by  expansion  (or  compression),  The  third  expresses  cooling  (or 
heating)  provided  by  the  space  gradient  of  the  electrons  "carrying 
away"  velocity,  The  fourth  member  expresses  cooling  caused  by 
heating  of  cold  electrons  (just  created)  to  a given  temperature 
(=incroase  of  their  kinetic  energy  to  the  mean  kinetic  energy  of 
the  remaining  electrons),  the  fifth  member  expresses  cooling  or 
heating  Mfft  of  electrons  under  the  influence  of  collisions 
and  /...illegible,../ 

a)cooting  by  elastic  collisions  with  heavy  particles,  b)  cooling 
by  charging  (or  stepwise  charging)  of  neutral  (or  charged,  or 
Ionized)  particles,  d)heatlng  by  collisions  of  other  kinds. 

T<et  us  write  furtheri 
(15)  A **  -N. pVJotf.  - 


where  cx  is  the  frequency  of  the  k-th  order  of  collision  for 

K 
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one  electron  and  q U,c  is  the  energy  consumed  (or  gained,  then  U,c 
will  be  negative)  during  one  such  collision  of  the  Ic-th  order.  Thus 
we  shall  obtain  the  equation  (1.4)  in  the  following  form 
(16) 


KI  /do  vO\  K,  _ , dN.O  n mYioU 

N.  (t~  + U'T'I  = - 5ir,  •— i N,G  - O — - N.  Z.A'JoL <.  • 

\et  dz  J dz  dz  6f  ‘ 

3. LINEARIZATION  OF  THE  ELECTRON  ENERGY  TRANSFER 

EQUATION  AND  DISCUSSION 

Let  us  assume  such  a plasma  disturbance  where  the  magnitudes 
in  the  positive  column  are  modified  in  such  a manner  that  their 


values  become  i N,  = Nt0  + nr , |n,|  « n,0  , 

0 = €>0  + 9,  |9|  < 0O, 

E = E0  4 e,  jej  = |£0|  , 


(17) 


da* 

a»  = »10  + -—.9 
SO 


(the  sub-index  o denotes  the  vaLuea  in  the  homogeneous  positive 

column).  Equation  (l6)  for  the  homogeneous  column  has  the  following 
form  i 
(18) 


9o«,N,0£0  - A',o  - 0 . 

A 


Let  us  substitute  the  relationships  (l?)  in  equation  (l6),  substract 
<%)  equation  (l8)  and  after  rearranging  we  obtaini 

^•T7  4 — + N,o  Y,  ~ = QoUr^'rOe  + <Jo".E0H,  + a„  a..  RW 

* cz  ‘ - j./.0  p - ",  IcWJoCL  - 0 . 

V in  ^ & 8' 

From  here  it  is  evident  thaty'f’he  analogous  equation  (l)  only 

the  second  and  third  left  side  members  and  the  first  member  on  the 

right  side  of  that  equation  are  talcen  into  account,  that  is,  the 

derivation  with  respect  to  time  and  the  members  related  to  the 

variable  count  of  particles  have  been  negLected. 

For  further  discussion  we  must  express  the  "carrying  away" 


velocity.  We  shalL  use  for  this  purpose  the  continuity  equation 
( 1 3 ) • Thus  it  is  valldi 


(where  D and  M,  are  diffusion  coefficient  and  the  rnobilitv  of 
e / e J 

electrons,  respectively). 


also 

(20) 


b,  tX  „ r 
- N.Vz  ' "-E  ■ 


Then  from  equation  (19)  it  follows  that  if  we  consider  De  as  a 
constant  and  disregard  various  furtfcrlz  further  causes  contribu- 
ting to  the  "carrying  away"  velocity  up  which  are  present  in 
equation  (20),  we  wilL  not  get  in  equation  (19)  any  member  expre- 
ssing the  ^diffusion  ^heat  of  electrons  caused  by  the  temperature 

ia  ( , 

gradient  (that  is,  a member  containing  c)  & f\  ) but  solely 
a member  expressing  the  electrons  thermal  diffusion  caused  by  the 
gradient  of  particles'  count  (that  is,  a member  containing 
a ^ &Z.  ).  The  member  containing  ° /Oz.  may  be  obtained 

/ s 

only  by  precising  the  expression  (20)  for  the  "carrying  away"  veloci- 
ty. T,et  us  consider  if  it  is  possible  to  maKe  this  equation  more 
accurate  by  these  two  methods i 

1.  By  expanding  up  by  one  more  member,  by  the  component  of  the 
''carrying  away"  velocity  related  to  the  electrons'  thermal  gradient 

(of  the  type  (-l/TG)  grad  D T0).  The  mechanism 
of  heat  diffusion  of  the  electrons  is  of  two  types  /4/  i 
a)  transfer  of  icinetic  energy  by  collisions  of  the  electron-electron 
type i but  because  these  collisions  taKe  place  in  a low  pressure 
plasma  at  small  currents  they  are  unLlKely  to  occur  and  may  be 
neglected,  b)  Penetration  of  electrons  from  higher  Levels  of  mean 
Kinetic  energy  into  lower  mean  Kinetic  energy  Levelst  this  manner 
of  diffusion  is  given  onLy  by  the  concentration  gradient  which  is 
included  just  in  the  "carrying  away"  velocity. 


-fe- 


2.  By  considering  the  dependence  on  temperature  of  De  (thermaL 
diffusion  - nee  //$/)  where  the  mobiLity^M.  e can  be  considered 
In  first  approximation  ar.  constant  (this  may  easily  be  generalized), 

(_/<-  ) ljr  a Linear 


also  ?>e  in  the  expression  i)e  =.  ^ 


/ 


function  of  temperature,  then 

1 dD.N, 


a.  = - /',£  . 

A',  dz 


thus 


du 

dz 


and  because 


i = Bi  p'hV  _ PJK  _ I dO,  dn.  Dc  d2n,  de 
■ 1*1  \to)  dz2  N;  dz  dz  A '.  dz2  ’ 

D-  - 


we  have 

du. 

^ & 


Dj  (du, V cD,  £** 

A?W  cr2 


J_  ?«,  _ 0,  _ a? 

A,  <?0  £>z  A/,  cz2  ^ ' dz 


From  this  we  see  that  the  fourth  member  on  the  right  side 
of  equation  (19)*  when  a Linear  dependence  on  temperature  is  con- 
sidered for  , indicates  besides  other  processes  also  the  electron 
thermaL  diffusion  under  the  influence  of  a temperature  gradient 
(second  member)  and  the  Dart  ides’  count  gradient  (fourth  member 

of  expression  (2i).  By  other  processes  we  mean  here  a dependence 

"temperature 

on  the  total  time  variation  of  the  electrons4 ’ on  the  gra- 
dient of  electric  fieLd  deviation  from  the  equilibrium  state,  on 
the  particLes’  count  gradient  and  on  the  mixed  relationships  with 
the  particles’  count  gradient  and  the  electron  temperature 
gradient. 

4.  CONCLUSION 

For  smalL  discharge  currents  (up  to^lOmA/cm2 ) we  performed 
in  a uni-dimensionaL  approximation  an  energy  averaging  of  the 
Boltzmann-Vlasov  icinetlc  equation  assuming  a quasl-Maxwe L lian 
eLoctron  velocities  distribution.  V/e  deducted  equations  which 
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correspond  with  equation  (i  ) and  have  considered  a method  for 
a further  approximation.  In  conclusion  the  author  than'cs  M.T.ichov 
C.Fc.,  V.Venely  and  T,,  ^e'carlcov  G.f.c.  for  their  valuable  suggestions 
and  V.Fuchsov  for,  besides,  a careful  proofreading  and  comments. 


ADDENDUM 

>leulatl»n  of  integrals. 

as  a consequence  of  independence  of 


1 * f c1  dC:  ‘ 

J(C)  Ct 


time,  coordinates  and  instantaneous  velocity  we  may,  after  solving 
for  re  (3)  and  (4)  and  substituting 
( ) Cj  = C cos  3 , c2  = C sin  3 cos  <p  , c3  - C sin  3 sin  ip 


write 


w he  re 


( D2  ) 


f c>  * dc  = d Nfy^y2  r r p 

et  dt  \ * / Jc-oJi-oJ^o' 

• c c*‘«*  cJ  sin  3 dC  d3  d<j> . e 

/i  = 2y,infi/,C. 


From  the  relationship  |i|  ^ < Wf>  (,btain 

the  magnitude  of  the  "carrying  away"  velocity  of  the  electrons j 
for  the  density  of  the  current  i ^ lOmA/cm2  and  concentration 
V -v,  1010  cm"3  we  obtain  u_  l06cm/s.  Further  for  T_^2-  lO1*  °K  we 

arrive  at  X s2/cm2.  From  here  we  see  that  the  factor 

Ug  o 


(D3) 


1 ! _ + ...il 

1! 


c 

By  carrying  the  integration  through c^and  substituting  cos O =n 
we  obtain 

f c2  Q'dC  2 it  A/,  ( ^-'V  ^ P C4c-’*"*c,f  e*’dsdC. 

J(f,  \ * / Jco  J.  — i 
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fi  2 v>  h 

e1"  ds  = sinh  li  - 2 X 
h 


oo  L 2 I 

X ■ 

j=o  (2/  1 1)! 


and  then  with  ( D2  ) ’ * f c2  <Yf  dC  = * 2 £ ..2'f 

J (^)  it  Ct  L 71  / ' °(2/  + ])  Jo 


C!luC""'C’dC, 


f c‘e!?6C.  tN,i2'^(2y,m,y  • „» 

Cl  it  1-0  (2/)!! 


a^ter  integration 

J(r,  ^ »*»  (2»)! 

solving  the  argument  for  the  L-th  member 

1 iiji  * (io15)'  1 .(io'1)2'  « io"3,+ 15 , 

that  in,  each  following  member  in  by  3 ordern  nmalLer  twan  the 
preceding  and  it  nufficen  to  Limit  ournelven  to  1=0.  Ir  we  intro- 
duce (6)  we  get* 

f e >cfr  dC  = -(  2 Nto). 

J<o  ct  ci\m,  J 


2,  2.  I c\c\:r)f,dC  i calcuLated  by  determining  f from  (3) 

J (O 

and  (4)  and  writing  out  the  ncalar  effectn  and  then  nubntitute 
in  (Dl)»  then  (D2)  and  (D3)  obtaining 


(7s  c-v.".c*  e*‘o*»  sjn  g cos  3 dCdfl d</>  -f 


0 

+ — 
0JC 


’2* 


C5C  y*",'c,cfceo'a  sin2  3 cos  dC  d3  d<p  4- 


ojojo 
0 r*>  r*  r2* 


+ | ||  C-e  ■'"'•cVt04<  sin2  3 sin  </>  dC  d.9  dip  > . 

dy  J o J o J o J 


Both  the  Latter  integration  membern  through  will  be  elimi- 
nated and  the  firnt  member  we  solve  fully  an  in  the  nrecedinp 
cane,  Thun  we  have 
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we  first  simplify  the  integration 


f f?(F'  V')s. 

I.  J(r>  \»'c  J 


ilC 


by  parts  (see  reference  /3/  for  detaiLn) 

J C’Cf  Vr)'/,dC--.lf  cF./.dC; 

J(r)  V».  / n,c  J (<•) 


talcing  rrom  the  re Lationshio  (9)  and  assuming  that 
( ^c.  — [c  x H]^  = o,  we  obtain  by 


out  the  scaLar  effects,  again  f*r  three  members,  of  which 
the  first  one  is  of  the  zero  tyre  (by  integration  through 
Thus  we  solve  again  as  in  the  case  l.  and  obtain 


-q0NtueE. 

m. 


RF.FF.RKNGF.r, 


|l)  Pi  kxrfk  L.,  Kfiiifi  V.:  C/cch  J Phys.  H I?  vl')62),  296. 
(21  Gu.wgv^ki)  V.  Flcklnieskij  lok  v y.wc  i. 

[3]  Kumii:  J.:  IJsod  do  tcorie  phzmatu  (uccbni  text). 

Hi  I liens  V.,  Visfi  v V.:  xoukromc  sdi’lcni. 
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